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Abstract 

'p ,■ It is shown that the magnetohydrodynamic equilibrium states of an axisym- 

c*2 ' metric toroidal plasma with finite resistivity and flows parallel to the magnetic 

field are governed by a second-order partial differential equation for the poloidal 
magnetic flux function ip coupled with a Bernoulli type equation for the plasma 
density (which are identical in form to the corresponding ideal MHD equilibrium 
equations) along with the relation A*ip = V c a. (Here, A* is the Grad-Schluter- 
Shafranov operator, a is the conductivity and V c is the constant toroidal-loop volt- 
age divided by 2ir). In particular, for incompressible flows the above mentioned 
partial differential equation becomes elliptic and decouples from the Bernoulli 
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equation [H. Tasso and G. N. Throumoulopoulos, Phys. Plasmas 5, 2378 (1998)]. 
For a conductivity of the form a = a(R,ip) (R is the distance from the axis of 
y* symmetry) several classes of analytic equilibria with incompressible flows can be 

constructed having qualitatively plausible a profiles, i.e. profiles with a taking a 
maximum value close to the magnetic axis and a minimum value on the plasma 
surface. For o = cr(ip) consideration of the relation A*ip = V c a(ip) in the vicinity 
of the magnetic axis leads therein to a proof of the non-existence of either com- 
pressible or incompressible equilibria. This result can be extended to the more 
general case of non-parallel flows lying within the magnetic surfaces. 
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I. Introduction 

In addition to the case of the long living astrophysical plasmas, understanding 
the equilibrium properties of resistive fusion plasmas is important, particularly 
in view of the next step devices which will possibly demand pulse lengths of the 
order of 10 3 sees (or more for an ITER size machine) (H and Refs. cited therein). 
Theoretically, however, it was proved by Tasso that resistive equilibria with 
a = ff(if}) are not compatible with the Grad-Schliiter-Shafranov equation. (Here, 
a is the conductivity and ip is the poloidal magnetic flux function.) The non- 
existence of static axisymmetric resistive equilibria with a uniform conductivity 
was also suggested recently || |J ||. Also, in the collisional regime Pfirsch and 
Schliiter showed that the toroidal curvature gives rise to an enhanced diffusion, 
which is related to the conductivity parallel to the magnetic field. In the above 
mentioned studies the inertial-force flow term p(v • V)v is neglected in the equa- 
tion of momentum conservation. For ion flow velocities of the order of 100 Km/m, 
which have been observed in neutral-beam-heating experiments |6|, [7|, || the term 
p(v ■ V)v can not be considered negligible. Therefore, it is worthwhile to inves- 
tigate the nonlinear resistive equilibrium, in particular to address the following 
issues: (a) the impact of the non-linear flow in the Pfirsch-Schluter diffusion, and 
(b) the existence of resistive equilibria, in particular equilibria with a = <j(i(j). 
Since the magnetohydrodynamic (MHD) equilibrium with arbitrary flows and fi- 
nite conductivity is a very difficult problem, in a recent study || we considered 
an axisymmetric toroidal plasma with purely toroidal flow including the term 
p(v • V)v in the momentum-conservation equation. It was shown that the non- 
linear flow does not affect the static-equilibrium situation, i.e a = cr(ip) equilibria 
are not possible. 

A way of constructing more plausible equilibria from the physical point of 
view could be by considering flows less restricted in direction. Taking also into 
account the fact that the poloidal flow in the edge region of magnetic-confinement 
systems plays a role in the transition from the low-confinement mode to the high- 
confinement mode, in the present report we extend our previous studies to the case 
of flows having non-vanishing poloidal components in addition to toroidal ones. 
Because of the difficulty of the problem we consider flows parallel to the magnetic 
field. Some of the conclusions, however, can be extended to non-parallel flows 
lying within the magnetic surfaces. It is also noted that possible equilibria with 



parallel flows would be free of Pfirsch-Schliiter diffusion because the convective 
term v x B in the Ohm's low vanishes. The main conclusion is that for the 
system under consideration the existence of equilibria depends crucially on the 
spatial dependence of conductivity. The report is organized as follows. The 
equilibrium equations for an axisymmetric toroidal resistive plasma with parallel 
flows surrounded by a conductor are derived in Sec. II. The existence of solutions 
is then examined in Sec. Ill for the cases a = a(R,ip) (R is the distance from 
the axis of symmetry), and o = <j(ip). Sec. IV summarizes our conclusions. 

II. Equilibrium equations 

The MHD equilibrium states of a plasma with scalar conductivity are governed 
by the following set of equations, written in standard notations and convenient 
units: 

V ■ (pv) = 0, (1) 

p (vV)v=jxB-VP, (2) 

V x E = 0, (3) 

VxB=j, (4) 

V-B = 0, (5) 

E + vxB = -. (6) 

a 

It is pointed out that, unlike to the usual procedure followed in equilibrium studies 

with flow [[10], 11], [12, [13|, [TJ|, [15] in the present work an equation of state is not 



included in the above set of equations from the outset and therefore the equation 
of state independent Eqs. ( IToD and ( |T5D below are first derived. This alternative 



procedure is convenient because the equilibrium problem is then further reduced 
for specific cases associated with several equations of state. 

The system under consideration is a toroidal axisymmetric magnetically con- 
fined plasma, which is surrounded by a conductor (see Fig. 1 of Ref. 0). With 
the use of cylindrical coordinates R, <p, z the position of the surface of the con- 
ductor is specified by some boundary curve in the (R, z) plane. The equilibrium 
quantities do not depend on the azimuthal coordinate (p. Consequently, the di- 
vergence free magnetic field B and current density j can be expressed, with the 



aid of Ampere's low (|j), in terms of the stream functions ip{R, z) and I(R, z) as 

B = /V0 + V0 x W, (7) 

and 

j = A*V>V0 - V0 x VI. (8) 

Here, A* is the elliptic operator defined by A* = i? 2 V • (V/-R 2 ) and constant ip 
surfaces are magnetic surfaces. Also, it is assumed that the plasma elements flow 
solely along B: 

pv = KB, (9) 

where K is a function of R and z. Acting the divergence operator on Eq. (|9|) and 
taking into account Eq. ([I]) one obtains VK • B = 0. Therefore, the function K 
is a surface quantity: 

K = K(iJ;). (10) 

Another surface quantity is identified from the toroidal component of the mo- 
mentum conservation equation (0): 

(l-y)/ = X(V0. (11) 

From Eq. flTT| ) it follows that, unlike the case in static equilibria, I is not (in 
general) a surface quantity. Furthermore, expressing the time independent electric 
field by 

E = -V$ + V c V(j), (12) 

where V c is the constant toroidal-loop voltage divided by 2n, the poloidal and 
toroidal components of Ohm's law (|6]), respectively, yield 

v* = Y12UL (13) 

a 

and 

A*<0 = V c a = E^Ro. (14) 

Here, E^ is the toroidal component of E. Eq. ( |i~4"|) has an impact on the boundary 
conditions, i.e. the component of E tangential to the plasma-conductor interface 
does not vanish. Therefore, the container can not be considered perfectly con- 
ducting. Accordingly, Ohm's law with finite conductivity applied in the vicinity 



of the plasma-conductor interface does not permit the existence of a surface layer 



of current [16]. It is now assumed that the position of the conductor is such that 
its surface coincides with the outermost of the closed magnetic surfaces. Thus, 
the condition B ■ n = 0, where n is the outward unit vector normal to the plasma 
surface, holds in the plasma-conductor interface and therefore the pressure P 
must vanish on the boundary. It is noticed that this is possible only in equilib- 
rium, because in the framework of resistive MHD time dependent equations, the 
magnetic flux is not conserved. With the aid of equations (f7|)- (|ll]) the compo- 
nents of Eq. (§) along B and perpendicular to a magnetic surface are put in the 

respective forms 

r (K 2 B 2 \ VPl , N 
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Eq. ([TBI) has a singularity when 
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1. 



(17) 



On the basis of Eq. 



for pv and the definitions v\ = - — %-*- for the Alfven 
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velocity associated with the poloidal magnetic field and the Mach number 
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Eq. ( p!7|) can be written as M 2 = 1. 

Summarizing, the resistive MHD equilibrium of an axisymmetric toroidal 
plasma with parallel flow is governed by the set of Eqs. (fT^), ( |ToD and flT6"|). 
Owing to the direction of the flow parallel to B, Eqs. (|T5| ) and (|1^) do not con- 
tain the conductivity and are identical in form to the corresponding equations 
governing ideal equilibria. Therefore, on the one hand, several properties of the 
ideal equilibria, e.g. the Shafranov shift of the magnetic surfaces and the de- 
tachment of the isobaric surfaces from the magnetic surfaces (see the discussion 



following Eq. fl26|) in Sec IIC) remain valid. On the other hand, as will be shown 
in Sec. Ill, the conductivity a in Eq. ([14]) plays an important role on the existence 
of equilibria. 

To reduce further equations fll5|) and (0), the starting set of equations (|])-(§) 
must be supplemented by an equation of state, e.g. P = P(p,T), along with an 
equation determining the transport of internal energy. Such a rigorous treatment, 
however, makes the equilibrium problem very cumbersome. Alternatively, one 
can assume additional properties for the magnetic surfaces associated with either 
isentropic processes, or isothermal processes, or incompressible flows. These three 
cases are separately examined in the remainder of this section. 

A. Isentropic magnetic surfaces 

We consider a plasma with large but finite conductivity such that for times 
short compared with the diffusion time scale, the dissipative term rs j 2 /a can be 
neglected. This permits one to assume conservation of the entropy: v ■ VS = 0, 
which on account of Eq. @ leads to S = S(ip) (S is the specific entropy). 
It is noted that the case S = S(tp) was considered in investigations on ideal 
equilibria with arbitrary flows Jll], [12|] and purely toroidal flows fL7 , 18| , as well 



as on resistive equilibria with purely toroidal flows |§. In addition, the plasma 
is assumed to being a perfect gas whose internal energy density W is simply 
proportional to the temperature. Then, the equations for the thermodynamic 
potentials lead to ]17| 

P = A(S)p< (19) 

and 

w = Wy.pT-1 = E m (20) 

7-1 7 

Here, A = A(S) is an arbitrary function of S, H — W + P/p is the specific 
enthalpy and 7 is the ratio of specific heats. For simplicity and without loss 
of generality we choose the function A to be identical with S. Consequently, 
integration of Eq. ( |15|) yields 

« 2 B 2 , 7 „-_! 



2p 2 7 
Eq. ( |16D reduces then to 

K 2 \ W 



Sp^ 1 = Hty). (21) 
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p R? 



+ (v • B)K' + ^-X' + P H' - p^S' = 0, (22) 



where the prime denotes differentiation with respect to ip. Apart from a factor 
1/(7—1) in the last term of the right-hand side ([1/(7— l)]p 7 S" instead of p" f S') Eq. 
(|2"2| ) is identical in form with the corresponding ideal MHD equation obtained by 
Hameiri [0 (Eq. (7) therein). It should be noted that Eq. (22) remains regular 
for the case of isothermal plasmas (7 = 1) while Hameiri's result would make 
the equilibrium equation strangely singular. In particular, for S = S(ip) and 



T = const. Eq. ( |i"9|) leads to p = p(ifj) and consequently the incompressibility 
equation V • v = follows from Eq. ([!]). Incompressible flows, however, are 
described by Eq. fl27D below which is free of the above mentioned singularity. 

Unlike the case of static equilibria, Eq. (E^) is not always elliptic; there are 
three critical values of the poloidal-flow Mach-number M 2 at which the type 
of this equation changes, i.e. it becomes alternatively elliptic and hyperbolic 
|fiO| , |T2|j . The toroidal flow is not involved in these transitions because this is 
incompressible by axisymmetry and, therefore, does not relate to hyperbolicity 
(see also the discussion in the beginning of Sec. IIC). 

B. Isothermal magnetic surfaces 

Since for fusion plasmas the thermal conduction along B is expected to be fast 
in relation to the heat transport perpendicular to a magnetic surface, equilibria 
with isothermal magnetic surfaces are a reasonable approximation JT7], [TS|, |T^, ^D], 



21, B3]. In particular, the even simpler case of isothermal resistive equilibria has 



also been considered 123 



For T = T{jp) integration of Eq. ( jT5|) leads to 

K 2 B 2 
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+ ATlnp = #(</>), (23) 



where A is the proportionality constant in the ideal gas law P = XpT. Conse- 
quently, Eq. fll6|) reduces to 



V 



— 'g 



+ (v • B)K' + ^-X' + pH' - Ap(l - logp)T' = 0. (24) 



We remark that apart from the fact that the S terms have been replaced by T 
terms, Eqs. (|23|) and (p4 ) are identical with the respective Eqs. (^l|) and (|22|) . 



C. Incompressible flows 



The existence of hyperbolic regimes may be dangerous for plasma confinement 
because they are associated with shock waves which can cause equilibrium degra- 
dation. In this respect incompressible flows are of particular interest because, 
as is well known from gas dynamics, it is the compressibility that can give rise 
to shock waves; thus for incompressible flows the equilibrium equation becomes 
always elliptic. For V ■ v = it follows from Eqs. (Jl|) and @ that the density is 
a surface quantity 

P = PW, (25) 

consistent with the fact that in fusion experiments equilibrium density gradients 
parallel to B have not been observed. 

With the aid of Eq. (p5), integration of Eq. (|15|) yields an expression for the 



pressure: 

v 2 K 2 B 2 

P = PsW -Y = Ps ~-2p-- (26) 

We note here that, unlike in static equilibria, in the presence of flow magnetic 

surfaces in general do not coincide with isobaric surfaces because Eq. (|2|) implies 

that B • VP in general differs from zero. In this respect, the term P s {ip) is the 

static part of the pressure which does not vanish when v = 0. If it is now assumed 

that — — y^ 1 and Eq. (pq) is inserted into Eq. (|16D, the latter reduces to the 

elliptic differential equation 

(1 - M 2 )A> - 1 -{M 2 )'\ VVf + \ [y^Jp) + R2p s = °- ( 27 ) 

Eq. ([27]) is identical in form to the corresponding ideal equilibrium equation 



(Eq. (22) of Ref. p2|). It is also noted that special cases of incompressible 



ideal equilibria have been investigated in Refs. [Q and |25|. Unlike to the 



corresponding sets of compressible S = S(ip) equations ([21]) and (p2|), and T = 
T(t/>) equations ( |23|) and (|24|), Eq. ( ^7f) is decoupled from Eq. (p6|) . Once the 
solutions of Eq. (|27| ) are known, Eq. ( ]26| ) only determines the pressure. 

III. The existence of solutions in relation to the 
conductivity profile 

We shall show that the compatibility of Eq. (0) containing the conductivity 
a with the "ideal" equations ( |15|) and fll6f) depends crucially on the spatial de- 
pendence of a. In this respect the cases a = a(R, if)), and a = a (if)) are examined 
below. 



A. a = a(R,iJj) 
An explicit spatial dependence of a, in addition to that of if), is interesting 



because it makes the equilibrium problem well posed, i.e. in this case Eq. (|14|) 
can be decoupled from the other Eqs. ([TJ) and ([H]). A possible explicit spatial 
dependence of a can be justified by the following arguments: (a) Even in Spitzer 
conductivity, a = aT^ 2 , the quantity a has a (weak) spatial dependence and (b) 
cylindrically symmetric resistive a = crNJ) equilibria are possible [[J and therefore 
the non-existence of axisymmetric static toroidal a = oNJ) equilibria is related 
to the toroidicity involving through the scale factor |V0| = 1/R; this could also 
imply an explicit dependence of a on R. In addition, we may remark that the 
neoclassical conductivity depends on the aspect ratio A because the fraction of 



trapped particles relates to A (see |26j and Refs. cited therein). It should be 



noted, however, that a knowledge of the cr-profile in the various collisionality 
regimes of magnetic confinement has not been obtained to date. 

For us the main advantage in allowing a = cr(R, ijj) lies in the fact that Eq. 
(|T4] ) can then be considered as a formula determining the conductivity 

A*<0 . , 

a = -^ (28) 

provided ip is known. Also, the poloidal electric field can then be obtained by 
Eq. (0). 

To determine ip in the case of compressible flows with isentropic magnetic 
surfaces the set of Eqs. ([H]) and fl22[) , which are coupled through the density p, 
should be solved numerically under appropriate boundary conditions. This can 
be accomplished by the existing ideal MHD equilibrium codes 0, [14], [y|. The 
problem of compressible flows with isothermal magnetic surfaces [Eqs. ( fffi) and 
(|24|)1 can be solved in a similar way. 

For incompressible flows ip can be determined by Eq. (p7|) alone, which is 
amendable to several classes of analytic solutions. In particular, sheared- poloidal- 
flow equilibria associated with "radial" (poloidal) electric fields which play a role 
in the L-H transition can be constructed by means of the transformation [27, 28 

UN)) = [ [1 - M 2 N)') 1/2 ] dnj}, M 2 < 1, (29) 

Jo 



Under this transformation Eq. fl2"T|) reduces (after dividing by (1 — M 2 ) 1 / 2 ) to 



A*U + 



1 d 



X 2 



2 dU \ 1 - M 2 



+ iT 



dU 



(30) 



It is noted here that the requirement M 2 < 1 in transformation fl2"9|) implies that 



vi < vi 



where v s = ('yP/p) 1 ' 2 is the sound speed. This follows from Eqs. (|T8|) 



and (in Gaussian units) 



i) ^ 2 >W 



1. 



Since, according to experimental evidence in tokamaks [29], the (maximum) value 
of the ion poloidal velocity in the edge region during the L-H transition is of the 
order of 10 Km/sec and the ion temperature is of the order of 1 KeV, the scaling 
v p <C v s is satisfied in this region. Therefore, the restriction M 2 < 1 is of non- 



operational relevance. The simplest solution of Eq. 
const., X 2 = const, and P s oc ■0, is given by 



corresponding to M 2 



V> 






R r J \Rr 



(31) 



where ip c is the ip value on the magnetic axis located at (z — 0, R — R c ) and d 



is a parameter related to the shape of flux surfaces. Equation (|31|) describes the 
Hill's vortex configuration [ 30| . The conductivity then follows from Eq. (|28|) : 



a 



Or 



R 
Rr 



RrJ \Rr 



(32) 



where a c is the value of a on the magnetic axis. The conductivity profile in 
the middle-plane z = is illustrated in Fig. 1. We remark the outward displace- 
ment of the maximum-conductivity position -R max with respect to R c (R max /R c = 
2/y / 3) and the asymetry of the inner part of the profile as compared with the 
outer part due to the explicit R dependence of a. 

B. 



a 



*M 



For this case we consider Eq. (|14[) in the vicinity of the magnetic axis by trans- 
forming the coordinates from (R,z,<f)) to (x,y,<f)) (Fig. 2). The transformation 
is given by 



R = R c + x = R c + rcosO 
z = y = —rsinO. 



(33) 
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The quantities ip(x,y) and a(ip) are then expanded to second-order in x and y: 

9 1 

x y 

i/>(x, if}) = ipc + ci— + c 2 — + c 3 xy + ... (34) 

and 

x y 

a = a c + o-i (^ -tp c ) + ... = cr c + °"i( c iy + c 2y + C 3 X ^ + •••) + •••• (35) 

Here, Ci = (d 2 ip/dx 2 ) ci c 2 = (d 2 %jj/dy 2 ) c , c 3 = (d 2 ijj/dxdy) c , a c is the conductivity 
on the magnetic axis and 0\ = const. On the basis of Eqs. ([34]) and (|35|) Eq. 
A*-0 = V c a(ijj) becomes a polynomial in x and y which should vanish identically. 
This requirement leads to C\ = c 3 = and, therefore, it follows from Eq. (34) 
that the magnetic surfaces in the vicinity of the magnetic axis are not closed 
surfaces. 

The non-existence of <j(i/j) equilibria with closed magnetic surfaces can be 
extended to the case of non-parallel flows lying within the magnetic surfaces. 
Indeed, if the relation v • V^ = is assumed instead of v || B, the toroidal 
component of Eq. (|6|) leads again to Eq. ( jTjj ) . 

A possible proof of the non-existence of rj = r)(ip) equilibria far from the 
magnetic axis has not been obtained to date. It may be noted, however, that 
for a = cr(ip), Eq. ( |16D becomes parabolic. This follows by considering in this 
equation the determinant T> of the symmetric matrix of coefficients. On account 
of A*-^ = V c a(ip), and p = p(R,ip, \Vip\) by Eq. (|i"5l), the second derivatives of 
equation flTE|) are contained only in the term 



T^W vm ' ■ ^ 



which comes from the term V ■ [(1 — K 2 / 'p)V^/-R 2 ]. Subsequent evaluation of 
T> leads to T> — 0. Therefore, the function ip is (over) restricted everywhere to 
satisfy a parabolic equation and the elliptic equation A*-^ = V c a(ip). 

IV. Conclusions 

The equilibrium of an axisymmetric plasma with flow parallel to the magnetic field 
has been investigated within the framework of the resistive magnetohydrodynamic 
(MHD) theory. For the system under consideration the equilibrium equations 
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reduce to a set of a second-order differential equation for the poloidal magnetic 
flux function ip coupled through the density with an algebraic Bernoulli equation, 
which are identical in form with the corresponding ideal MHD equations, and the 
equation A*ip = V c a. (A*, V c and a are the Grad-Schliiter-Shafranov elliptic 
operator, the constant toroidal loop voltage and the conductivity, respectively. 
The existence of solutions of the above mentioned set of equations is sensitive to 
the spatial dependence of a. 

For a conductivity of the form a = cr(R, ip), Eq. A*ip = V c a can be considered 
uncoupled to the other two equations, thus determining only the conductivity. 
For compressible flows and isentopic magnetic surfaces the differential equation 
for ip [(Eq. (|22|) 1, pending on the value of the poloidal flow, can be either elliptic 
or hyperbolic. Solutions of the set of this equation and the coupled Bernoulli 
equation [Eq. (pl|)1 can be obtained numerically. The problem of compressible 
equilibria with isothermal magnetic surfaces [Eqs. ( ^3|) and (|24]) 1 can be solved 
in a similar way. For incompressible equilibria ip obeys an elliptic differential 
equation [(Eq. (p7])1, uncoupled to the associated Bernoulli equation [Eq. (|26|)1 
which just determines the pressure. Several classes of analytic equilibria with 
incompressible flows having qualitatively plausible a profiles, i.e, profiles with a 
taking a maximum value close to the magnetic axis and a minimum value on the 
plasma surface, can be constructed. In particular, sheared-poloidal-flow equilibria 
can be derived by means of the transformation fl2"5|) for ip. 

For a = a(ip) appreciation of A*tp = V c a in the vicinity of the magnetic 
axis proves therein, irrespective of plasma compressibility, the non-existence of 
closed magnetic surfaces. This result can be extended to the case of non-parallel 
flows lying within the magnetic surfaces. In addition, for parallel flows ip is 
(over)restricted to satisfy throughout the plasma an elliptic and a parabolic dif- 
ferential equations. 

According to the results of the present investigation, the existence of resistive 
equilibria is sensitive to the spatial dependence of conductivity. Thus, the task of 
obtaining this dependence in the various confinement regimes of fusion plasmas 
may deserve further experimental and theoretical investigations. A conductivity 
with a spatial dependence in addition to that of if), on the one hand, would open 
up the possibility of the existence of several classes of resistive equilibria free of 
Pfirsch-Schluter diffusion. On the other hand, a strict Spitzer-like conductivity, 
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a = cr(ip), should imply the persistence of a Pfirsch-Schliiter-like diffusion also in 
the non-linear flow regime. 
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Figure captions 
FIG. 1. The conductivity profile on the middle-plane z = described by Eq. 

FIG. 2. The system of coordinates (x,y,(f>). 
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FIG. 1. The conductivity profile on the middle-plane z = described by Eq. (R2| 




FIG. 2. The system of coordinates (x, y, 



